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ABSTRACT 


The improved and completed mathematical model of the solar radiation force 
and torques, which has been previously developed for the Mariner 10 Venus/ 
Mercury spacecraft mission, is used for a detailed analysis of the effects of 
the solar light pressure on the Helios spacecraft. Due to the fact that the 
main body of the Helios spacecraft is a surface of enclosure, inside of which 
most of the reradiated thermal energy is lost, expressions for the portion 
of the solar radiation force, produced by the thermal reradiation, had to be 
given a different form. Hence the need for the derivation of a somewhat 
different theoretical model for the force acting on the main body of the space- 
craft . 
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INTRODUCTION 


I, 

The heliocentric orbit of the Helios spacecraft is an extended ellipse 
of eccentrictiy 0.54 . After the injection into its orbit at its aphelion 
point at approximately one astronomical unit from the Sun, the spacecraft is 
falling freely toward the Sun until it reaches its perihelion point at approxi- 
mately 0.3 astronomical units from the Sun. It moves in the orbital plane of 
the Earth (ecliptic plane) with an orbital period of approximately 180 days. 

The proximity of the spacecraft to the Sun in the region of its orbit 
near its perihelion point causes a relatively large effect of the solar 
radiation force on the reflecting surfaces of the spacecraft, thus creating a 
necessity for a detailed analysis of the effects of the solar pressure on the 
orbital motion of Helios. The mathematical and physical concepts of the model, 

used in this analysis, are essentially the same as those outlined in References 

1 through 6. 

The side viev/ of the Helios spacecraft is shown in Figure 1. The 
principal parts of the spacecraft are: 1) the main body of the spacecraft, 

consisting of two truncated cones, connected by a circular cylinder in the 
middle; 2) three antennae: high-, medium-, and low-gain antenna, and 3) 

four booms carrying instruments . More detailed side and top views of Helios are 

shown on Figures 2 and 3. 


*Numerical values used throughout this report are for Helios A spacecraft. 
JPL Technical Memorandum 33-778 



The mathematical model for the solar radiation force will be given in 
the form of a short outline. In addition to this, particular attention will be 
paid to the derivation of the expression for the thermal reradiation of the 
absorbed solar energy, emanated from the relatively large surface of enclosure 
of the main body of the spacecraft. 

Since the spacecraft is rotating at 60 rpm, the most suitable reference 
system, which will be used for the calculation of the solar radiation force, 
is an inertial ecliptic system, translated into the center of mass of the 
spacecraft. 


II. THE SOLAR RADIATION FORCE 

The solar radiation force acting on the spacecraft is the sum of solar 
radiation force acting on one particular component of the spacecraft, and is 
given by (Reference 1) 


F = - 


// 


[f(0)N + (1 - 3 y) U COS0] dS, 


( 1 ) 


where 

f(0) = 2 3 ycos^ 0 + B(f) [y(1 - 3) 4- (1 - y) K(r, 0) ] cos0 (2) 

“62 

Here = 4.51312 x 10 N/m is the solar constant, p is the heliocentric 
distance of the spacecraft, expressed in astronomical units, 0 is the local 
angle of incidence, N is the unit vector along the local normal to the irradiated 
surface, U is the unit vector along the spacecraft-Sun direction, B(f^ is the 
diffuse reflection coefficient, Y is the total reflection coefficient and 3 
is the specular reflection coefficient of the illuminated surface. If Lambert's 


2 
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cosine law of diffuse reflection is used 


B(f) = 


S is the irradiated surface area of the spacecraft's component, and 
K(r, 9) is the thermal reradiation function. From References 1 and 9, the 
explicit expression for K(r, 9) is 


K(r, 0) 


£ T £ 

F F B 


T + e 
F F B 



(3) 


in the case when both spacecraft's components are exposed to space, and T- 

F B 

are respectively absolute temperatures and e and e emissivities of the 

F B 

front (illuminated) and back surfacss of the spacecraft's component. From 
the balance between the total absorbed solar energy and the amount of energy 
reradiated from both surfaces of a piecewise slab structure, using Gauss' 
theorem and assuming that the areas of the front and back surfaces are the 
same, we have 




K (1 - r) 


a P' 


COS0 


(4) 


where a is the Stefan's constant 

a = 5.6697 x 10 ^ kg/sec^ °K^ , 

and 

K = 1.353 X 10^ kg/see^. 

For an elementary slab of thickness Laplace's equation and boundary conditions 
yield 
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( 5 ) 




° 4 

-£■ -f £_ T 

B k B 


where k is the thermal conductivity of the material. Equations (4) and (5) 
yield 


K(r^0) = 1 


n 

2kp" sec6 . . 

Kg 5,(1 - y) F B^ 


) 


( 6 ) 


this expression is linear with respect to the difference of absolute temperatures 
of the front and back sides. 

For surfaces of enclosure the absorbed thermal energy which is reradiated 
inside the body is multiply scattered in all directions, since the interior of 
the surface contains various parts of different shapes. The so-scattered 
particles lose their momentum and only an insignificant small portion is re- 
radiated outside. Thus, the expression for K(r,6) becomes 


K(r,0) = 


^F 




(7) 


where, as before, T„ and T have to be determined from Equations (4) and (5) 

r B 


Finally, for adiabatic surfaces. 


K(r,0) = 1, 


For surfaces for which T„ = T^ (no temperature gradient) , 

i? B 


K(r,0) = K = 

■F B 


( 8 ) 
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III. 


OF THE SPACECRAFT . 

The main body of the Helios spacecraft is a surface of enclosure inside of 
which the multiply scattered radiation is lost. Fol^wing the procedure in 
Reference 1 we can expand the function K(r,0) in the following manner. First 
we write the approximate expression for the temperature ratio 


AR(r,0) 

1 + BR(r,0) 
B 


1 + (A + B) R(r,0) 
1 + BR(r,e) 



Further expansion of the temperature ratio leads to the following 
expression 

■j ^ 2 

(^) = 1 + 4AR(r,0) +2A (3A - 2B) [R(r, 0) ] (10) 
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where 


3A - 2B = - 3A 


E — P 

^ 

£ + £ 

F B 


= - 3AK. 


Denoting by 


K, 


"p "b 


we obnain from Equation (7) the approximate expression 


K(r,e) = + PR(r,0) + [R(r,e)]‘ 


or, explicitly, 


« 3/4 o 3/2 

, . -,,/Cos0, , „.cos6^ 

K(r,0) - + Q( 

P P 


( 11 ) 


where 


<s + V' 


A, 


.3 * 

+ 63) 


( 12 ) 


IV. SOLAR RADIATION FORGE ON THE MAIN BODY OF THE SPACECRAFT 

The main body of Helios, whose side and top views are shora in Figure 4, 
consists of two truncated, sixteen-sided pyramids, connected in the middle by a 
sixteen-sided prism. On account of the very fast spinning rate of Helios, we 
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shall considerably simplify our calculations if we approximate the body of 
Helios by two truncated cones, connected in the middle by a circular cylinder. 
The average radii of the truncated cones are calcualted in such a manner that 
the areas, which play a major role in the expressions for the solar pressure 
effects, are preserved. In other words the radii are calculated so that the 
area of one truncated pyramid is equal to the area of one truncated cone, and 
the same rule is used for the cylinder in the middle. Thus, for the larger 
radius of the frustrum we obtain 

= 1. 372 m , 

and for the smaller radius which, at the same time is the radius of the middle 
cylinder, we find 


= 0 . 871 m 


(see Figure 5). The height of one of the frustrums is H - 78.65 cm, and the 
height of the middle cylinder is h = 55.0 cm. 


In order to facilitate the calculation of the solar radiation force 
acting on the main body of the spacecraft, we shall introduce a noniuertial, 
slowly rotating ecliptic reference frame ^1‘ The z^-axxs xs poxntxng 

toward the north pole of the ecliptic, x^ y^-plane is the ecliptic plane, and 
the y^-axis lies along the instantaneous spacecraf t-Sun direction. Incidentally, 
we should mention that, without any significant loss in accuracy, we can assume 


that the z^-axis is also the spin-axis of the spacecraft, 


In the X y z system the equation of the conic surface is 

111 


f(x^, yj_> 
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x^ ^ + y^^ - (z^ + 5)^ tan a , 


(13) 


7 



where 






h 

2 


is the z-coordinate of the vortex of the cone. Hence, the unit vector along 
the local normal to the surface is 


N = 


V(j) 
” V(J) 


cosa 

+ ? 


cota 
y. cota 

i 

(z^ + O tanay 


where a is the cone angle, given by 


tana = 


\ ~ ^2 
H 


or a = 32?75. 


Introducing cylindrical coordinates r^, A by 


we obtain 


^1 “ » 

= r^ sinA , 

^1 ~ ^^1 tana 


(cosX cosa 
N = I sinA cosa 
- sina 


(14) 
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Also, 



and 


ds = 


sina 


Since the main body of the spacecraft is perfectly symmetrical with respect 
to the three coordinate axes , the only existing component of the solar radiation 
force will act along the y^-axis . Integrating Equation (1) we obtain the y^- 
component of the force acting on both frustrums in the form 


21 H 
s 


YF 


3P 


- (R^ + Rj) {aBj, Y 


p cos^a +1 [Yp (1 - 3p) 


P (cosa) 
£ 


7/4 


+ (1 - Yp) Kp] cosa + 3(1 - Yp + U ~ Yp) ^1 

5/2 

T 

2 


(15) 


Q (cosa) 

+ (1 - I 


’F" .3 

P 


where Bj,, Yp. Qp. S SrustruTns, 

I^ and are given by 


n 

= y* (sinl)^^'^^ dl, I 2 = y (sinl)^^^ dX , 
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or, in terms of beta and gamma functions (Reference 10), 




= 2 


2 

11/4 r ^8 ^ 


r if-) 


= 1.382638 , 


^2 = !> 


2 9 

,7/2 ' V 


r (|) 


= 1.248598 


The complete integrated expression for the y^-compcnent of the solar radiation 
force acting on the surface of the middle circular cylinder can be taken from 
Reference 6. It is 


(16) 


where 6 , Y , P » Q > > ^.re the values of 6, Y, P, Q, and K for the middle 

c 'c c c c 

cylinder. The total solar radiation force acting on the whole main body of 
Helios is then 

p _ = p 4- ,p _ 

B YF YC ’ 

acting in the Sun-spacecraft direction. 

V. THE SOLAR RADIATION FORCE ACTING ON OTHER APPENDAGES OF THE HELIOS SPACECR^T 

A very rough calculation of the total solar radiation force acting on all 
other major components of the spacecraft shows that the maximum contribution to 
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+ 3(1 - B Y ) + (1 - T„) 
c c c 


'^1 ^ , "2 V 

, „3/2 3 I 


10 



the force, due to the irradiation of all those components, does not exceed 6.8 
per cent of the force exerted on the main body of the spacecraft. Still, in 
order to perform a complete analysis of the solar pressure effects on Helios, 
we shall furnish a short synopsis of the solar radiation forces acting on other 
parts of the spacecraft, shown on Figure 6 and explained in Table 1. 

All the components have simple geometric shapes for which the expressions 
for the solar radiation force have already been derived in References 1-6; 
they all can be considered as having purely adiabatic surfaces for which 
K(r,6) = 1. Particular attention has been given to the derivation of the solar 
radiation force acting on the despun high— gain antenna, pointed toward the 
Earth, and its partially shadowed frame. The force acting on the four rapidly 
rotating booms and instruments has been averaged over one period of rotation 
(1 sec.). It is very evident that due to the smallness of the forces involved, 
it has been very beneficial. 


COMPONENT NO. 

DESCRIPTION 

1 

Main body 

2 

High- gain antenna 

3 

Medium-gain antenna 

4 

Low-gain antenna 

5 

Instruments and the dummy 


counterbalance 

6 

Instrument carrying booms 

7 

Booms 


TABLE 1 

Principal parts of the spacecraft. 
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and convenient to use approximate expressions which considerably reduce the 
complexity of the mathematical apparatus applied and ultimately Implemented 
into the computer programme. 

VI. COMPUTER PROGRAMME FOR THE CALCULATION OF THE SOLAR RADIATION FORCE . 

The solar radiation force Is a function of the heliocentric distance of 
the spacecraft and, up to a lesser extent, of the position of the spacecraft 
relative to the position of the Earth. Thus, in order to compute a series of 
values of the solar radiation force over a certain Interval of time, we should 
first compute the positions of the spacecraft in its heliocentric orbxt. The 
relatively small magnitude of the force {the maximum total force is approx- 
imately 2.63 X 10“'* N or 26.3 dynes) has led to the conclusion that an 
unperturbed, elliptical orbit would yield a more than satisfactory accuracy for 

that purpose. The orbital elements, which served as input data for the 

calculation of this orbit, were obtained from the personnel of the German 
Geselschaft fur Welf raumf orschung by J. Peyn of the University of Hamburg. 

The computer programme, listed in Appendix 1, is entirely self-explanatory 
via its comment statements, so that additional descriptions and explanations 
would undoubtedly be superfluous. Although the use of only one space-fixed 
system of reference axes would be perfectly sufficient (such as the ecliptic 
system for 1950.0, initially used) , the programme also computes the components 
of the solar radiation force along the axes of the Earth-equatorial reference 
frame of 1950.0. It also gives the magnitude of the total force acting on the 
spacecraft. The forces are computed for each day (sero-hour GMT) during one 
complete orbital period of the spacecraft (approximately 190 days) , starting 
with the day of the insertion into the elliptic orbit, five days after the 

JPL Technical Memorandum 33-778 


12 


launch time; the output is given in Appendix 2. In addition to this, a plotting 
programme has been used to depict graphically the behaviour of the force 
components over a period of AOO days. The graphs are shown in Appendix 3. 
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TOP VIEW 



Fig, 6. Principal spacecraft components 
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The Computer Programme 
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ivCME.'iCLATUKL 

tPSF = EMISSIVITY CF THE ILUJMI.NATED SURFACE 
EP'^B = EMISSIVITY OF THE BACK SURFACE 

GA*^MA = REFLECTIVITY COEFFICIENT OF THE ILLUMINATED SU'RF-CE 
BETA = specular REFLECTIVITY PARAMETER OF THE SURFACE 
AIGMA = STEFAN-S constant (3,&6<;7E-C6 KG/SEC^<-*3.Df GK<i-*4) 
iOLAR ■= solar CCNSTANT (1.3b3E + 0S <G/SEC-it*3) 

P5I = the SUN-SPACECRA ft-earth angle IuPS) in RmDIANS 

DPS I = THE SUN-SPACECRAFT-EARTH ANGLE (EPS' IN DEGRESS 

SEP = SUN-EARTH-SPACECRAFT ANGLE IN RADIANS 

D-^EP =^- SUN-EARTri-SPACLCRAFT anGlE IN DEGREES 

,<1 = JA‘ACRSKI-S CONSTANT EPSP/ ( EPSF + EkSB ) 

alpha = CONE ANGLE OF THE MAIN BODY OF THE SPACECRAFT 

AX = THE SEMI-MAJOR AXIS OF THE SPACECRAFT-5 ELLIPTIC 

Axe = the^semi'-majcr axis of the elliptic orbit of the earth ( 

tCC = THE ORBITAL ECCENTRICITY OF THE SPACECRAFT 
ECCE = ECCENTRICITY CF THE EAR-^H-S ORE I T 
EpALN = OBLIQUITY OF THE ECLIPTIC 

INCL = INCLINATION CF THE ORBITAL PLANE TO THE ECLIPTIC 
PLANE CF 1950.0 

mode = NODAL ANGLE OF THE ORBITAL PLANE OF cjPACECRAFT 
OMEGA = ARGUMENT OF THE PER I APS IS OF THE SPACECRAFT-S ORBIT 
EOMEGA = ARGUMENT OF THE EARTH-S PER LAPS I:: 

THETA = TRUE ANOMALY OF THE SPACECRAFT 
ETHETA = TRUE ANOMALY CF THE EARTH 
M = mean anomaly of the SPACECRAFT 

ME = mean anomaly OF THE EARTH 

E = ECCENTRIC ANOMALY OF THE SPACECRAFT 

EE = ECCENTRIC ANOMALY OF THE EARTH 

MSTART = MEAN ANOKALY OF THE SPACECR*vFT AT THE TIME OF 
INITIALIZATION T=TSTART 

MESTRT = MEAN ANOMALY OF THE EARTH AT THE TIME OF INITIALIZATION 
T3TART = INITIAL TIME 

GM = GRAVITATIONAL CONSTANT OF THE SUN 

MEAN = mean orbital MOTION OF THE SPACECRAFT 

MEAME = .MEAN ORBITAL MOTION OF THE EARTH 

PER = ORBITAL PERIOD OF THE SPACECRAFT 

PERE = ORBITAL PERIOD CF THE EARTH 

N~ = NLM3ER OF POINTS ON THE TRAJECTORY 

TSTEP = TIME STEP (DAYS) 

RHO = HELIOCENTRIC POSITION VECTOR OF THE SPACECRAFT IN 
ASTRONCMICAL LNITS (AU) 

RHOE = HELIOCENTRIC POSITION VECTOk OF THE EARTH IN 
ASTRONOMICAL UNITS t AU > 

UR - SPACECRAFT-EARTH DISTANCE (KM) 

R1 = outer RADIUS or THE TRUNCATED CONES (MEtERS) 

P2 = RADIUS OF Th-E iMIDDLt: CYLINDER (METERS) 

ri = HEIGHT OF THE TRUNCATED CONES (METERS) 

HO = HEIGHT OF THE MIDDLE CYLINDER (METERS) 

= average RADIUS OF THE HlGH-s MEDIUM-. AND LOW-GAIN 

A.NTENNA DOOMS 

hA = total height of the THREE ANTENNA B-jOi- o 

THECRETICAL formulation FOR THE THER,mAL RE-RADIATIOK 

1. equation- 

EPSF*TFRONT**4 ETA'*^EPS8-s^TEACK*'*f4 = 

(SOLAR /SIGMA) *( (AU/R > ^'•^2 1 *( 1. -gamma ) *COS ( THETA ) 

(ENERGY balance, GAUSS" THEOREM) 


EE 

MSTART 

MiESTRT 

T3TART 

GM 

mtean 

MEAME 

PER 

PERE 

N 

TSTEP 

RHO 
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no n o o o n ri n n n ri rt nr) n n o n o n 


2. EGUATIGK- 

TFROKT = T8ACK + ( D*SIOKA#EPiB/COND) *TBACfC**A 
(BOUf^GARY conditions, lAFLACE-S equation) 

LpOCri = TSTART = 197A, DECEhSER Tj, 00 HRS, 00 MIN, OC SEC, 

X,Y,Z ARE ECLIPTIC COORDINATES OF TriE SPACECRAFT 
Xu,Y»,iw ARE EQUATORIAL COORDIMAtES OF THE SPACECRAFT 
XE,YE,2E=C ARE ECLIPTIC COORDINATES OF THE EARTH 

references- 

1. R.M.GEGRGEVIC, TLChNiCAL MEMORANDUM 33-49A, OCTOBER 1,1971 

2. K.M.GECRGEVlC, TECHNICAL MEMORANDUM 391-429, MARCH 30,1973 

3. .R.M.GEGRGEVIC, TECHiUCAL MEMORANDUM 33-698, DECEMBER 1,1974 

4. R.M.GEURGEVIC, PROGRAM -MVMSOL-, MAY, 1974 


SPEC I F I CAT I ONS- 


C 


REAL INCL,NODE 
REAL IXX,IYY,IZ2 

real MST ART, MZERO.MESTRT, MEAN, MEANE, MASS 


REAL i 1,12,141 

REAL M.{40G) ,ME(4uu) ,NCRM(28) 

DIMENSION THEtA(4UC) .EtHETA ( 4Gt' ) ,X(AC0) ,Y(’Ju),R(A0o) ,Z(A* 
<• xE(400) ,YE(AcC) ,RE(4C0 ) 

dimension ux(Aco) ,LY (AGO) , UK Ago) , - , 

dimension PSl (ACO) ,DPSl (AQO) ,RH0(AQ0) ,RH0E(400) ,TIME(400) 
dimension E( ago ) ,EE(40G) ,DTMETA(400) .DETHTA(4C0) 

UlMENSICN OM(AoO) ,DME(400) ,DE(40U' ,D££(4C0) 

DIMENSION UR ( 400 ) , FBX ( 400 ) , FBY (400 , FBZ ( 400 ) ,BACC ( 40G ) , 

* BFCRCE(4C0) 

D I MENS I. ON XQ ( 40 0 ) , YG ( 4 OO ) , ZQ ( 4C0 ) , SEP ( 4C0 ) »DSEP (400 ) , 
cFx ( 400 ). ,BFY ( 40^- ) jbFZ ( AOO ) 

DIMENSION FTX(AuO) ,FTY(A0G) ,FTZ(A00) 


0 ) 


DIMENSION FTMAG(4vjO) ,TACC( 400) 
dimension TFY1[4Gu) ,FAV(A0C) ,FLAT(4CC) 

DIMENSION AREA(28) 

DIMENSION SDIST(4GC) 

DIMENSION flag (28) 

DATA p 1/3. 141592654/ , DAY/ 864GC,/ ,SPE ED/2. 99792 5E+C6/ 

DATA 11/1.332638/ , 12/1.246598/ 

DATA Rl/ 1,373/ ,R2/0. 86713 / ,H/0,7865/ ,HC/C .550/ 

data DT/128. 2152777778/ 



NAME 1.1 ST /INPUT/EPsF , EPS B,GAMAN, SIGMA, solar, beta ,C0ND,ECCE,T start , 

1 T.STEP , AXE jEOMEGA , EPSLN ,GM ,MASS , AL',MZER0,N .DEPTH ,TLNCH , IXX » I YY » 

2 l"zz , I NCL , NODE, Ax , OMEGA , ECC ,M£STRT , NORM .AREA ,GAMSA ,GAM.CC , RA , HA , 

3 ABETA 
C 

RAO = 18C.0/PI 
70U0 READ (5, INPUT) 
wRITE(6,3CC0 ) 

WRITE (6, INPUT I 
V.RITE (6 ,'lOOC ) 

3UOO FORMAT ( 1H1.4X//// ) 

100 0 FOkiTA T ( 1 ri 1 , 1 X , / ) 

DO Sou J=1 ,2 8 

AREA(J) = AREA( J ) /ICOOO. 

NORM ( U ) = NORi'i ( J ) / R AO 
SOb CONTINUE 

ALPHA = ATAN( (Rl-R2)/H) 

MEAN - SGRT(GM/AX*^^3 ) ■ 

MEANE = SQRTiGM/AXE^f-K-S) 
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CMLAN = f'.ilANfDAY^^RAD 
DMilAf.E = MEANt<tDAY«-RA-: 

I F { MESTRT . GT . 360. ) MECTRT=MEoTRT-36C . 0 
PER = 2,C*PI/ (MEAM*DA':‘) 

PERE = 2,0*PI/ (MEANE*DAY) 

ETA = SGRT( ( l.'J + ECC)/' l.C-ECC) ) 

ETAL = SGRT( I l.G+ECCE;/( l.C-ECCE) J 
IF(CMLGA.LT.v;.C)ON;EGA = OREGA + 36C.C 
FGTAaT = MZERO - DFEAC'!'--DT 
IFtMSTAKT.LT.O. InSTAR '=MSTART+36C . 

0 o 3 R I T £ ( 6 » 2 ) 

/. R I r E ( 6 » 3 ) AX » ECC , OMEGA ,HSTAR T , DMhAN , PER » I NCL » RODE 
ARITE(6»A) 

V‘,RITE(6.31)AXE,ECCE,EuWEGA,MESTRT,DMEANE.PERE 

2 FCRIIAT (36X. -HELIOCENTRIC EGLATORIAL ORBITAL PARAMETERS OF THE -» 

* -SPACECRAFT-//) 

3 FORMAT (38X, -AX =- , E16 . 8 , 1 X ,-KM-/ » 3?X »-ECC =- , F 1 6 . 1 1 / . 3 5X . -OMEG A =- 

1 .F16.11,1X*-DEG-/,3AX,-MSTART =- . F 1 6 . 1 1 » 1 X .-DEG- / , 36X . -ME AR 

2 F16. 11. 1X.-DEG/DAY-/.3AX. -PERIOD =- . F16 . 1 1 . IX .-DAYS-/ . 36X . 

3 -INCL =-.F16.11 .1X.-DEG-/.36X.-NODE =- . F 1 6 . 1 1 . IX . -DEG-/ / / ) 

A FORMAT! /36X .-HEL J OCEMR I C ECLIPTIC ORdITAL PARAMETLRS OF THE -. 

1 -EARTH-//) 

3i FORMAT (38X, -AX =- . E 16 . H . 1 X , -KM-/ . 3 7X . -ECC =- , F 1 6 • 1 1 / . 3 5 X . -OMEGA =- 

1 . F16.il, lx. -DLG-/.3-^X.-MSTART , F 1 6 . 11 , 1 X .-DEG-/ , 36X , -MEAN 

2 F16.il. IX, -D£G/DAY-/,34 X, -period =- , FI 6 . 1 1 , IX , -DAYS- / / / / ) 

OMEGA = OMEGA/ RAD 

E'OMEGA = EOMECA/RAD 
EPSLN = EPSLN/RAD 
SE = SIN(EPSLN) 

CE = CCS (EPSLN) 

MSTART = MSTART/RAD 
MEGTKT = MESTRT/RAD 
I NCL = I NCL/ RAD 
NODE = NODE/RAD 
Cl = COS (I NCL) 

SI = SIM INCL) 

CN = COS (NODE) 

SN = SIN (NODE) 

CO = COS (OMEGA) 

SO = SIN (OMEGA) 

CONVERSION OF EuUATORIAL ORBITAL PARAMETERS FOR 1950.0 INTO ECLIPTIC 
PARAMETERS FOR 1950. 

SC = CE*SI - S£#CI*CN 
SS = SE*SN 

SNI = SQKT(SS**2+SC**2) 

CSI = CI*CE + SI^^SE^CN 
INCL = ATAN2 (SNI , CSI ) 

DELCM = ATAN(SS/SC) 

OMEGA = CMEGA-DELCM 
SNODE = . SI*SIN( DELOM) /SE 
NODE =ASIN( SNODE) 

CINCL = INCL*RAD 
DOMEGA = OMEGA-K-RAD 
DNODE = NODE*RAD 
WRITE (6, 97.) 

iVRlTE(6.:98)DOMEGA,DINCL.DNODE 

97 F0RMAT(36X, -ecliptic ANGULAR CpaiTAL PARAMETERS OF THE SPACCC-, 

* -RAFT-//) ■ 

9S F0RMAT(35X, -OMEGA =- , F 1 6 . 1 1 , 1 X , -DEG- / , 36x I MC l =- , FI 6 . 1 1 , 1 X , 
-DEG-/ ,36X , -MODE =- , FI 6 . 1 1 , Ix ,-OEG-/ / /.) 

COMPUTATION OF UNPERTURBED POSITIONS OF THE SPACECRAFT AND THE EARTH 

31 = SIN (INCL) 

Cl - COS (I NCL) 
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3K = bINtMODE) 

CM = COS (NODE) 

SO = SLN( OMEGA) 

CO = COS (OMEGA) 

PX = CN*CO - S.\*SO*Cl 
PY = SN^^CO + CN*50*CI 
PI = SO*SI 

QX = -CN^^SO - SN*CC<fCI 
GY = -SN^SO + CN^^■CO*CI 
Q2 = CO*SI 
RX = SN*SI 
RY = -CN*SI 
KZ = Cl 


TESTl 

= 

PX**2 

+ 

PY*^f2 

+ 

PZ^Hf2 

-1 

TEST 2 

= 

QX*42 

+ 

GY**2 

+ 

QZ**2 

-1 

TESTS 

= 

RX«-J<-2 

+ 

RY^*-»2 

+ 

RZ^^*2 

-1 

TESTA 

= 

PX**2 

+ 

QX«-»2 

4* 

RX<^^2 

-1 

TESTS 

= 

PY**2 

+ 

OY^^*2 

+ 

RY<:-^2 

-1 

TESTS 


PZit*2 

+ 

QZ^<-*2 

+ 

RZ**2 

-1 

TEST? 

= 

PX*QX 

+ 

PY*GY 

+ 

PZ«-OZ 


TESTS 

= 

PX^<-RX 

+ 

py«-ry 

+ 

PZ*RZ 


TEST9 


GX*RX 

+ 

QY*RY 

+ 

GZ*RZ 



TESTIC = PX^<-PY + OX-^GY + RX^*RY 

TESTll = PX^<-PZ + QX*OZ + RX*RZ 

TESTiZ = PY»PZ + QY*G2 + RY*RZ 

WRITE (6. 8) 

rtRITE(6*909) PX.PY»PZ.QX»QY,0Z,RX.RY,RZ 
V> R ITE(6,11) 

WR1TE(6,8C12 )TEST1 ,TESTZ,TEbT3,TEbT4,TEST5,TEST6,TEST7,TEST8, 

* TEST9 ,TEST10, TESTll, TEST12 

8 FORMaT( /35X >-VECT0RS P»Q,R OF THt SPACECRAFT-// ,36X ,-ECL IPT IC-/ ) 
9u9 format (10X,-PX =-,FlC.6,l0X,'-PY , FI 0 • 6 , lOX ,-PZ =-,Fl0.6/, 

1 luX,-GX , Flo • 6 , lOX , -GY , F 10 • 6 » 1 CX » -GZ , FI 0« 6/ » lOX , 

2 -RX =-,F1C»6,10X,-RY =-,El0.6,l0X,-RZ =-,F10.6//) 

11 Format (//36x,-ORTHOGOMALITY TESTS-/) 

8012 F0kMAT(/(4F15.6) ) 

MP = N/SO + 1 
DO 16 NPT5 = I ,N 

M(ImPTS) = MSTART + (NPTS-1 I^MEAN^tTSTEP 
ML(iNPTS) = MESTRT + ( NPTS-1 ) *-MEANE*TSTEP 
EtXPTS) = AN0M(ECC,M(NPTS) ) 

EE(NPTS) = ANCM(ECCE,ME(NPTS) ) 

R(KPTS) = AX*(1.0-ECC*C0S(E(NPT5) M . 

RE(NPTS) = AXC*-(1.0-ECCE*COb( EE(NPTS) ) ) 

RHC(NPTS) = R(NPTS)/AU 
RHOE(NPTS) = RE(i\PTS)/AU 

THETA(NPTS) = 2.C*ATAN(ETA*TAN(E(NPTS)/2.0) ) 

ETEIEtA ( KPTS ) = 2. 0*AtAN ( ETAE*TAN ( £E( MPTS ) /2 . 0 ) ) 

1F( E(RPTS) .Gb.Z.O^^-pn E(NPTS ) = E(NPrs)-2.G*Pi 
IF (Ec(inPTS) .GE.2.u*Pl ) EE ( NPTS ) =EE ( NPTS ) -2 .0*P I 
IF(H(NPTS) .GE.2,0*Pl )M(NPTS) =M(NPTS)-2.0^PI 
I F (ME ( HPTS) . GE. 2. U*P I >ME ( NPTb )=ME{ NPTS ) -2 . 0-»P I 
XORb = R(MPTS)*C0S(THETA(NPTS) ) 

YORB = R(NPTS)*SIN(THETA{NPTS) I 

X(KPTS) = PX*XORB + QX^YORB 
Y(NPTS) = PY*XORB + QY*YORB 
Z(RPTS) = PZ*X0R3 + QZ*YORB 
XG(NPTS) = X(KPTS) 

YC(NPTS) = Y(NPTS)»CE - Z(NPTS)*SE 
ZG(KPTS) = Y(NPTS)*SE + Z(NPTS)*CE 

XE(NPTS) = RE(NPTS)*COS(lTHETA(NPTS) + EOMEGA) 

YE(NPTo) = Rn.(NPTS)*SlN( lTHETA(NPTS) + EOMEGA) 
time (NPTS) = TSTART + ( ( NPTS-1 ) *ToTEP ) /DA Y - TENCH 
UX(NPTS) = X(NPTS)-XE(NPTSf 
UY(NPTS) = Y(NPTS)-YE(NPTS) 

LiZ(NPTS) Z(NPTS) 
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URoQ = UX(NPT5)**2 + UY ( NPTii ) »«-2 + UZ(NPT5)*«-2 
UR(MPT5) = 5GRT(UR5Q) 

C 

C CX»CY*CZ ARE: COSPOiNhNTS OF THE CROSS-PRODUCT (VECTOR-PRODUCT) (R X R£ ) 

C 

C SiN(PSl) AM' SI.Y(SEP) ARE DETERMINED FROM CROSS-PRODUCTS (RE X R) AID 
C ( K X RE) , RESPECTIVELY. 

C CUS(PSl) AND COS(SEP) ARE DETERMINED FROM DOT-PRODUCTS (SCALAR PRODJZTS 
C (UR.R) AND (UR. RE), RESPECTIVELY. 

C 

C TnE SIGN OF SIN(PSI) IS DETERMINED IN SUCH A MANNER THAT SIN(P5D IS 
C POSITIVE IF CZ IS NEGATIVE. THE ANGLE PSI IS MEASURED FROM R TO UR IT 
C THE POSITIVE DIRECTICN 

CX = -Z(,NPTS ) *YE(NPrS) 

CY = ZUmPTS )*XE(NPTS ) 

CZ = XINPTS )^<-YE(NPTS) - Y ( NPTS ) *XE( NPTS ) 

CSQ = CX*^^2 + CY**2 ■’r CZ**2 

DOT - bX ( NPTS )*X( NPTS) + U Y ( NPTS ) *Y ( NPT S ) + UZ ( NPTS ) *Z ( NPTS ) 

spsi = s(;rt(Cso)/(ur(npts)*P(Npts) ) 

I F(CZ.LT .0. ) SPSI=-SP31 
CPSI = DOT/(UR(NPTS)*R(NPTS) ) 

PSI (NPTS) = ATAN2 ( SPSI ,CPSI ) 

DPSKNPTS) = PSI (Np-'-M «-RAD 

IF(DPSI ( NPTS 1 .LT.O ,.PS I( NPT S ) =DPS I ( NPTS ) +3 60 . 

PDOT = UX ( NPTi) )*XE ( NPTS) + i- Y ( NPTS ) *YE ( NPTS ) 

SMSEP = SPSI *R( NPTS )/RE( NPTS) 

COS£P.= -PDOT/(UR(NPT£)*RE(NPTS) ) 

SEP (NPTS) = ATAN2(SNSEP,C0SEP) 

DSEP(NPTS) = SEP ( NPTS )*RAD 

CALCULATION OF THE SHORTEST DISTANCE FROM THE SUN TO THE GEOCENTRIC 
LIGHT RAY PATH 

DIST = R(NPTS)*SIN(PSI (NPTS) ) 

SDIST(NPTS) = ABS(DIST) 

DO 18 LIST=1,NP 
LL = 50^t(LIST-l)+l 
IF (NPTS. EG.LL ) WRI TE( 6,180 ) 

130 F0KMAT(1H1»1X» 

1 -TIME-,6X,-X(i<K)-,9X»-Y(KM)-,9X,-Z(l<M)-,9x»~RH0(Ai„‘)-, 

2 5X,-XE(KMi)-,8X,-YE(<iM)-,10X, -RHOE- , 8X , -PS I -/ , 

3 IX, -(DAYS)-, b9X,-(AU)-,7X,-(DEG)-/ ) 

18 CONTINUE 
WRITE (6, 19) TIME (NPTS) , 

1 X(NPTS) ,Y(NPTS) ,Z(NPTS) ,RHO(NPTS> ,XE ( NPTS >, yE ( NPTS ) , 

2 RHOE(NPTS) ,DPSI (NPTS) 

19 F0RMAT(1X»F6.2,3E14.6,F12.6,2E1A.6,F12.6,F11.2) 

16 CONTINUE 

DO 33 NPT5=1,N 

DTHLTA(NPTS) = THETA(NPT5)*RAD 
OETHTA(NPTS) - ETHETA(NPT5)*RAD 

IF(DTHETA(NPTS) .LT.0.0 )DTHETA(NPTS)=DTHETA(NPTS)+360.0 
IF(DETHTA(NPTS) .LT.u.O)DETHI A(MPT:j)=DETHTA(NPTS)+260.0 
DM (NPTS) = ,M(NPTS )-*«-RAD 
DM£(NPT5) = HE(NPTS)*RA0 
DE(NPTS) = E(NPTS)<^RAD 
OEE(NPTS) = EE(NPTS)*RAO 

IF(DM(NPT5) .LT.O, 0)DM( NPTS) =DM( NPTS) +36 C. 

IF(DME(NPTS) .LT.O. 0)D.v,E( NPTS) =DME( NPTS) +360.0 
IF(DE(NPTS) .LT.O. 0) DEI NPTS) -DE( NPTS) +360. 

1F(DEE (NPTS) .LT.O. 0)6lE( NPTS )=DEE( NPTS 1+360.0 
IF(DE (NPTS) .GE.36(J.0)DL( NPTS) =DE (NPTS) -360.0 
IF(CEE(NPTS) .GE.360. 01 DEE ( NPTS) =DEE(, NPTS) -360.0 
IF(Dy.(NPT5).GE.360.C)0;'('IPTSi-DM(HrT£)-360.0 
IF(DME(NPT5) ,GE.3 60.0)DMe( NPTu)=:DMC( NPTS)-360„0 
DO 119 LIST=1»NP 
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LL = 50*(LI5T-1)+1 

IF(NPT5.tQ.LL)WRITE(6,lii2) 

18 2 Format ( 1H1,29x»-T I me-, 5X* -theta-, 6x,-C-*8X,-M- »5X »-TH1iTA^-, 5X, 

1 -EE-,7x,-ME-,7x,-5DIST-/,?9X,-(DAYS)-,6(4x,-(DEG)-) ,6x,-(KM)-/) 

WR^TE^t's^l^O) T IME( NPTS) , DTHEtA ( MPTS ) ,DE { NPT5> ,DM ( NPTS > , 

1 DETHTA ( KPTS) ,2>EE ( NPT5 > ,DME ( MPTS) tSDlST ( NPTS ) 

190 F0RMAT(29X,F6.2 ,5F9.2,E14.6 ) 

33 CONTINUE 

AUXILIARY parameters FOR THE CALCULATION OF THERMAL RE-RADIATION EFFECTS 

ELES = -lOCuOUU .o*-solar/speed 
EPS = EPSF+EPSB 
<1 = EPSF/EPS 

-SA- MEANS SOLAR ARRAY (TRUNCATED CONES) 

-CC- I-iEaNS central COMHARTMENT (CYLINDER) 


C 

C 

C 

r 

r 

c 

c 

c 

c 

c 

c 

c 


temsa = solar*(1.-gamsa)/(sigmay;-eps) 

TEMCC = S0LAR*( 1 .-GAMCC) / ( S IGMA^EPS) 

TSTSA=TEMSA^f*0.25 
TSTCC=TEMCC**0.25 
C = SIGMA*DEPTH*EPSB/CQND 
ASA=C*TGTSA**3. 

ACC = C*TSTCC**3 
SSA=3.*ASA*i<l 
BCC = 3.*ACC*K1 

PSA = 4.*EPSF»EPSd*ASA/EPSF^^*2 
PCC = 4.*EPSF*L?SB*ACCAEPSF**2 
QNUM = -2.0*EPSF*EPSE*(11.0^^EPSF-3.Q*EPSB) 

GSA = QNUM*( ASA-fr*2 ) /EPS**3 
QCC = 0NUM*( ACC4<-*2 ) /EPS*’" 3 

the body OE the HELIOS SPACECRAFT CONSISTS OF- 

1. MAIN BODY OF THE SPACECRAFT - TWO TRONCATED CONES AND THE CYLINDER IN 
THE MIDDLE. ONLY THIS PART CONTAINS THE THERMAL RE-RADIATION EFFECTS. 

2. THREE ANTENNA BOOMS- HIGH-, MEDIUM— , . AND LOV.'-GA IN ANTENNA BOOMS. 

3. FOUR ROTATING SPACECRAFT BOOMS. TnE SOLAR FORCE WILL BE CALCULATED 
BY AVLKA()ING OVER ONE PERIOD OF ROTATION. 

4. ANTENNA WIRES. 

3. FLAT PARTS OF THE FRAME OF THE HIGH-GAIN ANTENNA REFLECTOR. 

£. PARABOLICALLY SHAPED FRAME OF THE HIGH-GAIN ANTENNA. 

calculation OF THE SOLAR RADIATION FORCE ACTING ON THE MAIN bCOY CF^THE 
HELIOS SPACECRAFT 

1. CALCULATION OF THE SOLAR RAC 1 AT I ON FORCE ON T.-.O TRUNCATED CONES 


C 

C 


CTl = 2.'*^ELES*H*(R1+R2)/3. 

CA = COS (ALPHA) 

Tl = 4,*bETA*GAMSA*CA**2 

T2 = , 5*P 1 *CA* ( CAMS A* ( 1 . -BET A ) + ( 1 . -G/TMSA ) *tCl 
T3 = 3.*(1.-EETA*GAM5A) 

T4 = ( 1.-GA,MSA)*PSA*I1*CA**1.75 

T5 = ( 1.-GAMSA)*QSA*12*CA**2.5 

TRMl = GT1*{ Tl + T2 + T3) 

TRM2 = CT1*TA 
TRM3 = CT1*T5 

2. SOLAR RADIATION FCRCLON THE MIDDLE CYLINDER 
CT2 = 2.*£LES*R2*HC/3. 

j5 = 4 ,*bETA*MAMCC + 3 «* ( 1 . -^ETA*GAMCC ) 

T7 = . p*PI* ( GAMCC* ( 1 .-beta ) + ( 1 ."GAMCC )*‘(C1 ) 

T8 = (l.-GAMCC)*PCC*Il 
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c 

c 


c 

c 


T9 = ( l.-CAMCC)*GCC*12 
TRKiA = CT2^(T6 + T7) 
TRM5 = CTZ^^TS 
TRM6 = CT2^^T9 


TOTAL HcLlCO eCDY RADIATION rC.’XL 
TRUNCATED CONES PLUS THt FORCE ON 


— E vj i~ R pI A L' I /“. T I o * « F o 1^ C E 
THE NIDJLE CYLINDER 


ON 


THE TWO 


TERYil = TRHl + TRMA 
TERY.2 = TRM2 + TRM5 
TERMS = TRM3 + TRM6 


CALCULATION CF THE SOLAR RADIATION FORCE ON THE ANTENNA uCOM 


CTA = 2.*£LES*RAR'HHA/3. 

TIL = 4 • ■jpA3ETA*SAMCC + 3 •■*'■( 1 • ~A3nTA'**''jAi'iC'- ) 

Til = ,5*?I* (SAMCC^^ ( l.-AEETA} + ( l.-CAMCC )-ii-Kl ) 

TRM7 = CTA*( Tlu+Tll ) 

C CALCULATION OF THE SOLAR RADIATION FORCE ACTING ON THE FOUR ROTATING 
C spacecraft looms » AVERAGED OVER ONE PERIOD 

r 

ABl = .043 
AB.2 = .0254 
AB3 = .GC335 
AEHl =3.768 
A5H2 =1.58 
ADH3 = 14.0 
BG = .60 

AECF = ELES^fP I* ( 1 .-BG)-*!-! 2 ,*ABl-»i-ABHl + 2.«AP2*A2H2 + AE3*A5H3)/6. 

C CALCULATION OF THE SOLAR RADIATION FORCE ACTING ON THE ANTENNA REFLECTOR 
C WIRES 


WCF = -.2 
C 

DALPHA = ALPHA*RAD 
ELAM = ELES/lOQOOOO.O 
ATR2 = TERMI/1000000. 

ATR35 = TERM2/1000000. 

ATR5 = TERM3/10CCGOO. 

WRITEtO.lOOlDALPHA.ELESjELAM.Rl.TSTSA.TSTCC.PSA.PCC.QSA.QCC, 

ATR2.ATR35.ATR5 

WRITE( 6 . 6 ) 

luo format( 1 Hi,////,iix>-phy 5 Ical Constants for calculation of the - 

1 -SOLAR RADIATION FORCE ON THE f/lA I M-/ , 1 2X , -BODY OF THE HElIOS - 

2 -SPACECRAFT-. 

3 /////, 24X.-ALPHA >-,F12.6»1X»-DEG-/>25X»-ELES=-» 

4 E15.6/ »23X .-LAMBDA =-.E15.6/. 

5 27X .-X1 =-.F12 . 6 / .24X ,-T5TAR =- , 2 F9 . 3/ . ? 8 X . -P = - ♦ 2c 1 ^ . 5 / , 2ex 

6 -Q =-,2E14.5//,26X,-TP2 =-,E1A.5/,25X,-TP35 , E14 . 5 / , 2 6X . 

7 -TR5 =-.E14.5//) 

6 FORMAT (IHl. 25 (/) .25X. 

1 -COMPONENTS OF THE SOLAR RADIATION FORCE ACTING ON - 

2 -THE MAIN, BODY OF THE-/ »25x »-HEL IQS SPACECRAFT-. 

3 - IN THE SPACE-FIXED. BODY-CENTERED ECLIPTIC REFERENCE FRAME.-// 

4 .25X.- (FORCE -. 

5 -GIVEN IN E-06 NEWTONS. ACCELERATION IN E-11 K,M/SEC'!f*2 ) -/ / ) 
AVTOT = O.C 

DO 1 NPTS=1.N 

RMl = l.C/RrlO(NPTS) 

r- 

C TOTAL FORCE IN THE SUN-SPACECRAFT DIRECTION 

C' ■ ... 

C 1.- FORCE ON THE MAIN BODY OF THE SPACECRAFT 

FYl = TERMl*R.Ml->n<-2 + T£RM2*RMi**3 . 5 + TERM3*RMl^''-^5 
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r-i n n nri n f't n r>n n nr^i n n 


FOrtCt Cl\ TKh ANTCNNA □COM5 

PyA — y. ■» f._< r* 1 ^ O 

3.- F0i<CL ON THE-' FOUR kOJATING tJOOiMS OF ThE SPACECRAFT (AVERAGED CVFR 

ONE PERIOD) 

FYA5 = ABCF*R;'.1»*2 

A.- FORCE ON THE V.'IRES OF THE ANTENNa REFLECTOR 
FYA.v = V.CF«R.'mf*2 

total force (SUN OF FORCES !• THROUGH A.) IN THE SUN-SPACECRAFT DIRECTION 
FYTOT = FYI + FYA + FYAS + FYAW 

Fax, FBY, FBZ ARE CC.NPONENTS OF THE iOLAR RADIATION FORCE ACTING ON THE 
,.1/UN oODY OF the HELIOS SPACECRAFT, ALONG THE AXES OF THE SPACE-FIXED, 
^ODY-CENTERLD EClI-’TIC reference FRA'^IF 


Fex(NPTS) = 
FSYCNPTS) = 
FBZ(NPTS) = 


FY1*X(NPT0)/R(N'^TS) 

fyi^yinpts ;/r(np'^s) 

FY1'»‘Z(NPTS)/R(NPTS) 


GFX, bPY, and BFZ are COMPONENTS OF THE SOLAR RADIATION FORCE ACTii»o 
CN^THE MAIN BODY OF THE SPACECRAFT IN T-E INERT 3 Al EOUATCRIAL, ECDY-C^NTEo.^D 

reference frame 

DFX(NPTS) = FSX(NPTS) 

bFY(NPTS) = FBY(NPTS)*CE - FBZ ( NPTS I«-SE 
5FZ(NPTS) = FBY (NPTS)*SE + FEZ(NPTS)*CE 

BFORCE(NPTS) = ABS(FYl) 

BACC(NPTS) = 100.a*bFCRCE (NPTS) /M-'^’SS 
DO 5 LIST=1,NP 
LL = 50^^(LIST-1) +1 
IF(NPTS,E0.LL)WRITE(6,7) 

7 FORMAT ( IHl ,3C'X ,-T IME-, 9X , -FBX- , 1 1 X ,-FBY- , 1 IX , -FBZ- , 10 X , -F wR'-u. 

* 10X,-ACC-, 11X,-RH0-/) 

5 CONTINUl 

WRITE(6,9)TIME{NPT5) ,FBX(NPTS) ,FBY(NPTS) ,F3Z(MPTS) ,3FORCE(NPTS) ,, 

* 3ACC( NPTS ) »HHO(NPTS) 

9 FORMAT (29X,F6. 2 ,6F14.5 ) 

("•'I nil ATlOf'i 'F TH“ COMPONENTS OF THl SOLAR RAD I AT I ON FORCE ACTING 0?i THt. 
FRAME cJ ?he1nTENNA^^^ S-acE-FIXED, BODY-CENTERED REFERENCE 

frame 

: xtgt: =0. 

YTOT = C. 

ZTOT = C. 

do: 12 J = 1 ,28 

FLAGU) - C. 

U = OMEGA + ThLTA(NPTS) 

UPAP = U + NORM(J) + PS I (MPT S) 

ENX = -C05(UPAP) 

ENY = SIN (UPAP) 

PPA = P5I (MPTS) +. MORM( J) 

CPP = COS (PPA ) / 

CRF = ELES*RM1^<-*2*AREA( J)^*-CPP 
CFA = 2.*&LTA-*-GAI';AN-»CPP + 2.->:-(l.-aETA* 

IF ( DPS I ( NPTS ) .GT. 0 , .AND. DPS I ( N.PTS * .LT . 

IFIDPSI (NPTS ):.GT .U. .ANC.DP51 (NPTSH .LT. 

IF (DPS I (NPTS ) .GT. A5. .AND.DPS I ( NPTS) .LT 
IF (DP Si (NPTS ).GT. AS. .AND .OPS I ( NPTo ) .LT 
I F ( DPS 1 ( NPTS / .GT .AS . .AND .DPS 1 ( NPTis ) . LT 
j p ( ppij I ( NPTS ) • GT • ^0 • • AND . DPS I ( NPT S ) . LT 
IF (DPSI (NPTS ) •GT.90. .AND.DPSI (NPTS) .i_T 
I p ( DPS I ( NPTS ) .GT . 90 . . AND.DPSI ( NPTS ) . LT 
1 F ( DPS I ( NPTS ) .GT.90. .AND.DPSI (NPTS) . I,T 

IF ( DPSI ( NPTS ) .GT.9&. .AND.DPSI ( NPTS) . LT 


GAM AM )./3 . 
leCi )FLAG(1)=1. 
180. )FLAG{2)=1. 
.225.)FLAG(3)=1. 
.225. )FLAG(A)=1. 
.225.)FLAG(5)=1. 
.270. )FLAG(§ )=:. 
.225 . )FlAG( 7 )=1. 
.130. )FLAG(8)=l. 
. 18C. ) FLAG(9 )=1. 
.270. ) FLAG( 1C) =: 
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^ n ri n n nr\ nn n r\r. n n r>n 


1 F ( DPS I ( !\PTS ) .bT . 1 3b . . . oP ? I ( hPTS ) . L T . 3 1 5 . ) F L i G ( 1 1 ) = . 

IF(JPSI ( ,\PTS ) .uT. 3 35../'-^iu.DFSU-\PTS) .LT.?3b. )-LAC-( 12 . 

IF(DPoI ( i-'iPTS) .5T.13b..A'lD.DPSI ( T S ) . LT . 3 1 b . ) "L AG ( I 3 ) . 

1 F (DPs I ( FiPTS ) .CT . 1 35. . ANiJ.CPSI ( .’4PTS ) . LT . Z 70 , 1 =LAC-f lA ) = ' . 

I F I DPS I ( ■‘■iPTS ) .GT . 1 . C.P^ 1 ( PPTS ) . LT .2 70 . ) PlAG( 1 5 ' = . 

1F(DPSI ( I'iPTS) »GT.IBo..'‘'Nl.«uPsI {?';PTSJ.l-T .PTC. )-lAG ( 16)= » 
IF(DPS1 ( i.-PTS 1 .GT. 1 BG. . Ar;D.DPCI (hPTS) .L T .350. ) •-■lAC.( 17 ) = . 


1 F (DPo I ( iSPTS ) .GT . 18^. . . A-MD.''. - si (PPl J .LT .3 b( . ! 'LAG( 1 ? ' = 

IF (DPSI ! NPTs ) .sT. . AfiL- .PPG I (r.'PTS J 
1 .Ai'.D.l,P 5 I { KP7b ) .LT. 36C. ) FLA ; ( 19 ) =1 . 

IFIDPSI ( NPTS ) .5T.0. .AN;:, DPS I <FPTS) 

1 .AAD.uPG 1( ^i°TS) .LT.360. ) FLA 2 3 ) =1 . 

IF(DPSI ( NPTS ) .GT.G. .AND. DPSI (APTS' 

1 .At, D. DPS I ( MPTS) .LT. 36c. ) FL^ C : 21 ) =1. 

IF (DP5l ( f'PTS ) «GT »v;» .AtJD.DPbl (APTS' 
l.Ai.D.DPSl (APTS) .LT.360. )FLAG(22) = 1. 

I F (DpsI t ,'<PT? ) .C-T. C . • A.NP • 'DPs I (f’PTS A 
I.AmD.DPSI CnPTS) .LT.360. ) FI, A G ( 2 3 ) = 1 . 

IF(DPSI ( iNPTS ) .GT. 0. .MND.i-PS I (APTS) 

1. ADD. DPS I ( DPTS) .LT.360. )FL A.'(24) =1 

I F (Dps I { MPTC ) .Gt. . . AM .IPf I ( \'PT o' 

I.AMD.DPSI (MPTS) .LT.360. ) FLA :,(?'; )=1. 

IF(DPSI (MPTS ) .GT.C. .AND.DPSI 
I.AMD.DPSI (,MPTS) .LT.36C. )FLAC(l6l=i. 

IFtDPSl (MPTS ) .GT.u. .AND. DPSI (MPfo' 
l.Ai’iD.DPSI (MPTS) .LT .360. )FLAo(27)=l. 

I F ( DPS I ( MPTS ) .GT . 0 . .AMD. DPS I ( f;P rs ) 

I.AMD.DPSI (MPTS) . L T . 360 . ) FLA G ( 2 8 1 = 1 . 

AFX = CRF* { CFA^-EiNX - ( 1 . -BE TA^-CA.-'.AK ) 4<-X ( NF T£ ) /^ ( MPT 3 ) ) »Fl AG ( 

AFY = GRF* ( CFA3(EMY - { 1 . -3ETA*GAMA''| ) «-y ( MF TS ) / 3 ( MPT d )) / 

AFZ = -CRF*( 1 .-BETA#GA,MAM ) ((2 (MPTS' *FLAG („ ) /R (NPTS ) 


LT 

. AS . 

• o • 

LP,j: 1 ( NF-T 

) • G” 

.225 . 

LT 

• AS . 

• C • 

DFSI ( '-'PT 

) . G T 

.22 5 . 

LT 

.45 . 

.OR. 

DPS 1 ( MPT 

. ) .GT 


LT 

.54. 

. OR • 

DPSI (MPT 

) .GT 

.233. 

LT 

. 9q • 

.Ok. 

Dpc I (MPT, 

. ) . C-T 

.21''.. 

LT 

• ?c. 

• CR. 

DPS' (MPT. 

, ) .GT 

.270. 

LT 

. 1-1 i 

..O’ 

. Dt.'; I ( '';r> 

S 1 .0 


LT 

.135 

. .CR 

.DTC I ( Mp 

' ) .G 

-."'15 „ 

LT 

.13 5 


.rp' r ! .’.iP 

S ) .0 

-.315 .. 

LT 

.13 5 

. .OR 

.''PS I ( \P 

S) .0 

-.315 . 


XTOT = XTOT + AFX 

YTOT = YTOT + AFY 

ZTOT = ZTOT + AFZ 

12 COiNTIMUE 
AFTX = XTCT 

AFTY = YTOT 

AFTZ = ZTOT 

CAvLCuLATlOM OF THE SOLAR RAD lAj ICN FORCE ACTING OK TFE^ PA'ABCL^C AfMTENf-A 
FRAME 


CF = .3*R,Ml«-^i-2/SlN(Psr(NPTS) ) 

XPT = uX (NPTS)^«•C0S (PSl (NPTS ) )/;jR(\PT- ) -X(.NPTS) / R(f PTS) 
YPT = LY ( MPTS )*COS ( PSi ( NPTS > > /UR ( MPTS ) -Y ( MPTS > /P ( '■ PTS ) 
PFX = -CF*XPT*.25*SIN(PSI (M^TS) ) 

PFY = -CF((-YPT*C0S ( PSI (MPTS ) ) 


TFX > TFY » 
acting on 
ECLIPTIC . 


TFZ, a, re the CoyPOMLNTS DF Thrl RESULTANT SC LAP >AD I A - IC'’-. ' DPCE 
the' whole body of' THE HEuIOS SPACt CRAFT, ALOMG hS A ,ES OF ThP 
SPACE-FIXED, BODY-CENTERED REFERFf C£ FRAME. 


TFX = FYTQT*X(NPTS)/R(NPTS) AFTX+PFx 
TFY = FYTOT*Y(NPTS )/R(NPT5) AFTY + PFy 
TFZ = FYT0T*Z ( NPTS ) /R( NPTS ) AFl Z 


TTX(NP.TS), FTYtiNPTS), FTZ(NPTS), ARE THE COMPONENTS CF THE RESJLTAM” 

SOLAR RADIATION FORCE ACTING CM THb SHOLE BODY OF THE HELIOS S-^ACEC-RAFT, 

ALONG THE axes OF THE LCUATORIAL, -SHACE-F I Xi D»' EODY-C ENTERED R-'FFRENCE FRAVE. 


FTX(NPTS) = TFX 

FTY(MPTS) = TFYP^CE - TFZ*SF 

FTZ(MPTS) ='TFY*sE ^ TFZ*CC 


't h'c MAviNltoDt: CF the Rt SULTAHT SCLAR RAOIAtIOM FORCE, IN E-o6 '•)£■, •■TOMS , IS 
FthAG ( MPTS ) = SORT (f TXCNPTS )**2 + FT Y ( HP 1 S ) ^'-^2 F F TZ ( MP TS ) ^'■*2 ' 
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C THE ACCELERATION OF THIS FORCE IN E-H KV,/SEC#*2 I£ 

TACCINPTS) = IOC. *FTMAG( NPTS) /MASS 

TOTAL FORCE IN THE SPACtCRAF.T-S'JN DIRECTION 

TOOT = FTX(NPTbH^X«(NPTS) + Fiy ( NPtS > *YO( M^TS ) 

* + FTZ(NPTS)*ZGti\PT6) 

TFYl(NPTS) = -TDOT/R(NPT5) 

lateral FORCE IN THE DIRECTION PERPENDICULAR TO THE_ sun-SPACECrAFT 
DIRECTION 

XEG = XEINPTS) 

YEC = YE(NPT3)*CE 

ZEQ = YE(NPTS)*5E 

RRE = X(NPTS)*XE(NPTS) + Y ( NPTS ) * YE ( NPTS ) 

DNK = UR(NPTS)*SIN(PSI (NPTS ) )^fR(MPT5)**2 

XIX = (RRE^XU(NPTS) - XEw>“-R ( NPTS ) ) /DMM 

XlY = (RR£<i-YQ( NPTS ) - YEC^R ( NPTS ) **2 > /DNN 

XIZ = (RRE^^ZQINPTS) - ZEL»R t NPTS ) «->2 } /DNM 

FLAT(NPTS) = (XIX^^FTX(inPTS) + XI Y^^FT Y ( NPTS 1 + X 1Z*FTZ ( NPT 5 ) ) / 

1 FTMAG(NPTS) 

AVERAGE VALUE OF ThE FORCE AT 1.0 AU 

FAViNPTG) = TFVl(NPTS)*-RH0(NPTS)**-2 
AVIDT = AVTOT + FAV(NPTS) 


1 CONTINUE 

AVTOT = AVTOT/N 
v\'RITE(6»2C0) 

ZjO FpRMATtlHl.25(/ ) .25X.-CCMPONEMTS OF THE SOlAR RADIATION FORCE - 

1 -ACTING ON ThE MAIN BODY OF THE -/ ,25X.-HEL lOS SPACECRAFT IN - 

2 -THE SPACE-FIXED ECUATORIAL, BODY-CENTERED - 

3 -reference SYSTEM. -/j?5X»-F0RCE - 
A -GIVEN IN E-06 NE'.vTONS-// ) 

DO 201 NPTS=1.N 
DO 202 LIST=1.NP 
LL = 50*(LIST-1 )+l 
IF (NPTS.EG.LL) WRI TE(6.203 J 

20! 3 FORMAT ( IH 1 . 1 X , -TI ME- . 6X > -X ( XM ) - . 9X » - Y I XM ) - , PX . -Z ( KM ) - » 9X , 

1 -KHO(AU)-',8X,-XFORCF-.BX,-YFORCt-,BX,-ZFORCE--,8X,-SEP(D)-/ ) 

202 CONTINUE 

WRITE(6»20A) TIM£( NPTS) ,X0( NPTS) , YO ( MPTS ) , ZD ( NPTS 1 ,RI10( MPTS) » 

1 bFX(KPTS),BFY(NPTS),bFZ.<NPT3) ,DSEP(NPTS) 

2oA FORMAT (F7e2.3E14.6»F12.5»AFlA.6) 

201 CONTINUE 

205 ?ORWn;Hl!25(/),25X.-COMPOKENTS OF THE RESULTANT “LAR RADIATION 

1 FORCE-/ ,25X, -ACTING ON THE WHOLE BODY OF THE HFLlOS oPACE-. 

2 -CRAFT-/ .25X, -ALONG THE AXES OF THE EQUATORIAL, SPACE-FlXED.- 

- body-centered-/ » 25X»-REFERENCE FRAME, THE MAGNITUDE OF THE^-. 
A -RESULTANT FORCE ,-/ ,2bX , -AND ITS f ’ 

5 -IN E-06 newtons-/, 25X, -ACCELERATION IN E-11 KM/SEC-*^^2-/ ) 

bo 296 NPTS=1 »N 
DO 297 LIST=1,NP 
LL = 50* (LIST-1 1 +1 

IF(NPTS.EQ.LL)WRITE(6,298) 

298 FORMATdHl, -TIME-, 5X,-PSI (D)-,AX,,-RH0(AU)-,8X, ..(-p. g. 

1 -FXEQU-, 9X,-FyEQL-, 9x,-FZEQU-, 9X, -FORCE-, 9x , -ACCEL- , 9X , 

2 -FY1-»11X,-FLAT-/ ) 

WRlTE dd99) TIME(NPTS) , dps: (NPTS) »RHO(NPTS) , 

1 FTX(NPTS) ,ETY(NPTS) ,FTZ(NPTS) ,FTMAG(MPTS) ,TACC(NPTS) , 

2 TFYKMPTS) ,FLAT(NPTS) 

299 F0RMAT(F6.1 ,F10.2,F11.7,7F1A.5 ) ^ - 

296 CONTINUE ..d'.::* , . 

QUALiTr 
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STOP 

END 


FUNCTION ANCM(ECCiM) 

THIS FUNCTION SUbROUTINE SOLVES THE KEPLER-S EOUATICN By ITERATIONS 
C 

real ,■■1 

DATA EPS/.OUoOOS/ 

ANCM = K 

2 ANCM = M + ECC^fSIN(ANOM) 

TEST = ANUM - K. - ECC*S I N ( ANOM ) 

IFtABS (TEST) .GT*EPS)GG TO 2 
IF(AdS (TEST) . LE. EPS ) RETURN 
END 


-» INPUT 


EPSF = 

0.7 4 

CP Sd = 

0.5c 

0 A.M SA 

= .65 

GAMCC 

= .85 

GAM AN 

= .95 


BETA = .95 
AdETA = .85 
SIGMA = 5.6697E-G6 
SOLAR = 1.353E+03 
COND = 0.25 
DEPTH = 0.01 
HA -■ 2.095 
RA = .03125 


ECC = 

.521807390542 

ECCE = 

(3.01672 

T5TART 

= C.O 

TSTEP 

= 86400, 

AX = . 

96301973563E+08 

AXE = 

0.1495989194E+9 

EP 3 LN 

= 23.44578889 

OMEGA 

= 257.4449C8784 


ECMEuA = 102.50966 
INCL = 23.A469891399 
NODE = .061380A058997 
MESTRT = 3 AO. 7615 
MZEKO = 71 .6320A12103 
TLNCH = -5. 

AJ = .1493978930E+9 

GPi = .132712A99390802500E + 12 

MASS = 356.9 

N =40 

I XX = 343.8 

I YY = 199.8 

IZZ = 382.4 

NORM = 90.,90.,135. ,135. »135.,18C.,1S0. , 135.,90.,I8C.,2^5.,225.,223 . . 
18u.*270.»225.^2708*270.^315.j3l5..3T5.j323.tC*?0»>36..45.>45.>45. 

AR EA = 140 . » 2 2 . » 1 60 » » 1 5 6 . f 15 6 8 f 1 40 . . 6 5 . » 1 9 . » 20 t > 2 0 . » 16 0 » » 1 56 . » 1 56 8 ) 6 5 . » 

i4(./8»198t2G.»<i2.»16u8»1568.156.>‘-'5.«14v*8.38..95.8l60 8>13688l5A8 

$END 

-FIN NIF- 
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APPENDIX 2 


Components of the Solar Radiation Force in the 

Earth-Equatorial Coordinate System of 1950.0 and 

_6 

its magnitude, given in 10 Newtons. 
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I me 

fx 

FY 

rz 

force 

5.0 

-‘4.^7 125 

-22.94794 

-9.95092 

25,4q9q7 

6,0 

-‘4.179^3 

-23.03033 

-9.98679 

25,44716 

7.0 

-3.87685 

-23. 1 1 72 1 

- 1 0.0246 1 

25.49369 

8.0 

-3.57008 

-23.20062 

- 1 0 . 0 6 '1 4 0 

25,54867 

9.0 

-3.27778 

-23.30461 

-10.10617 

25.61217 

lO.O 

-2,'?7557 

-23.40525 

-10. 14996 

25,68426 

1 1 .0 

-2.67106 

-23.51060 

-10.19580 

25,76504 

1 2 . 0 

-2,36386 

-23.62075 

-10.24372 

25.65462 

1 3 , 0 

-2,05357 

-23.73577 

-10.29375 

25,95314 

i^i.Q 

-1.73976 

-23.85576 

-10.34594 

26,06075 

l5,0 

- \ .82201 

-23.98080 

- 10 • 4 Q 0 3 2 

26,17762 

l6.0 

- 1 . 0 7 9 0 6 

-24. 1 1099 

-10.45693 

26,30394 

l7.0 

- . 7 7 2 8 ‘4 

-24.24642 

-10.51583 

26,43991 

1 8.0 

- , 8 .’4 □ ^4 7 

-24.38720 

- 1 0 * 577q4 

26,58577 

l9.0 

-.10230 

-24.53339 

-10.64060 

26,74173 

20,0 

, 2 '4 2 3 5 

-24.68516 

-10.70659 

26,9q812 

2 1 .0 

,57396 

-24.84259 

-10.77504 

27,08521 

22*0 

,9531*4 

-25.00579 

-10.84598 

27,2/331 

23 , 0 

1 .32050 

-25.17486 

-10.91949 

27,47276 

2*^.0 

1.67671 

-25.34992 

- 1 0 .99559 

27,68393 

25.0 

2.0825*4 

-25.53112 

-11.07436 

27,9q728 

2^.0 

2.87857 

-25.7 1849 

-1 1 .15501 

28,14314 

27.0 

2.88568 

-25.91219 

- 1 1 .2400 1 

28,39202 

20 . 0 

3 , 3o*47o 

-26 . 1 1 234 

- 1 1 .32702 

28,65444 

2‘?.0 

3 , 7 3 6 7 *4 

-26.3 19q1 

- 1 1 . 'U 6 8 6 

28,93088 

30.0 

‘4,18 20*4 

-26.53240 

- 1 1 .50962 

29,22206 

3 1.0 

*4.6/4243 

-26.75258 

- 1 1 . 60533 

2’. 52057 

32.0 

5,1 1870 

-26.77966 

-11.70404 

29,85109 

33,0 

‘i , 6 1 2 11 

-27.21379 

-1 1 .80582 

30. 19046 

3*+.0 

6,12371 

-27.45497 

-11.91 066 

30,54733 

35 , 0 

6,65'-497 

-27.7Q336 

-12*01864 

30,92269 

36.0 

7 . 2o735 

-27.959q5 

-12.12979 

31.31751 

37,0 

7 .78232 

-28.22208 

-12.24414 

31,73278 

38.0 

8,38179 

-28.49262 

-12.36176 

32.16961 

39,0 

9.0q7/}7 

-28.77068 

- 1 2.48265 

32.62978 

‘40,0 

9,66133 

-29,05630 

- 1 2.60683 

3 3 , 1 1 4 0 7 

‘4 1.0 

10,34561 

-29. 34956 

-12.73434 

33.62427 

‘42,0 

1 1 ,06239 

-29.65030 

-12.86514 

34.16167 

*43,0 

1 1 ,8 1*130 

-29.95080 

-12.99925 

34,72078 


12.70038 

-3O.202o6 

-1 3. 1 3982 

35,36888 

*45,0 

13.53340 

-3D.6o974 

-13.28232 

36,Oo735 

*46,0 

14.41022 

-30.94511 

-13.42818 

36,66201 

*47,0 

15,33441 

-31 .20805 

-13.57734 

37,39559 

‘4 8,0 

16,3o933 

-31 .63823 

-13.72966 

30,15069 

‘49,0 

17.33924 

-3 1 .99542 

-13»88504 

38.95061 

50,0 

1 9 . 5 1 0 3 5 

-32.90291 

-14.27910 

40,83073 

B 1 ,0 

20.65375 

-33.33328 

-14*46628 

4 1 .79663 

52.0 

21.06110 

-33.77050 

-14.65646 

42.81553 

S3.0 

23.13800 

-34 .2 1 344 

-14.84915 

43.89104 

B'i .0 

24.48956 

-34.66048 

-15.04366 

45,02665 

55,0 

25,92231 

-35.10909 

-15.23923 

46,22666 

56,0 

77.44345 

-35.55957 

-15*43496 

47,49594 
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time 

fx 

FY 

FZ 

force 

ri7 ,0 

29.06025 

-36.00676 

-15.62966 

48,83924 

58,0 

30.70111 

-36.44829 

-15.82195 

50.26220 

59,0 

32.61525 

-36.88033 

-16.01018 

51,77102 

6 0.0 

34,57281 

-37.29831 

-16. 1 9236 

53,37261 

6 1*0 

36.66475 

-37.69670 

-16.36612 

55.07445 

62, Q 

3 a . 9 0 3 9 2 

-38.Q69o4 

-16.52665 

56,88553 

6 3.0 

M 1 . 30292 

-38.4073? 

- 1 6 .67650 

5 8 , 8 1 4 6 5 

6^,0 

43.97747 

-38.70243 

-16.80567 

60.87316 

65,0 

46,64296 

-36.94292 

-16.91125 

63,07223 

66,0 

49.61669 

-39.11539 

-16.98744 

b5, 42479 

6 7,0 

52.01920 

-39 .20382 

- 1 7.02720 

67,94578 

6®.0 

56.26734 

-37 , 18886 

-17.02240 

70,65085 

69,0 

59 .9B636 

-39.04736 

-16,96276 

73.55810 

70,0 

6 3 . 9 9 8 1 /| 

-38.75136 

- 1 6 ,83625 

76,68695 

71.0 

60.29344 

-38 . 1 1 689 

-16.56309 

79,94515 

72,0 

72,95040 

-37,56014 

-16.32727 

83,67122 

73,0 

77,98418 

-36.57659 

- 1 5 .89957 

87 , 59096 

78,0 

6 3 , 3 9 1 n a 

-35.24900 

-15.32634 

9 1 ,82302 

75,0 

39,19863 

-33.51372 

- 1 4.57642 

96.39521 

76,0 

95,41034 

-3 1 . 28482 

-13.61265 

101.33461 

77,0 

102,05564 

-26.45895 

-12.39017 

106,67138 

78 , 0 

109,10184 

-24,91362 

- 1 0 . 8 5 5' 8 S 

1 12,43553 

79,0 

1 16,52925 

-20 .50495 

-8.94729 

1 18,65738 

80,0 

124,28747 

-15.06347 

-6.59096 

1 25,37035 

8 1,0 

132.2877 1 

-8 .39760 

-3.70373 

1 32 ,6q57 1 

8 2,0 

140,39789 

-,29066 

-.19157 

140,39832 

83,0 

148.42263 

9,48770 

4.04556 

1 48 ,78058 

8'4,0 

156,09394 

21.17564 

9.11 097 

157.78699 

85 , 0 

165.19121 

39.72110 
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APPENDIX 3 


Plots of the Components of the Solar Radiation Force 

in the Earth-Equatorial Coordinate System of 1950.0, 

—6 

and its magnitude, given in 10 Newtons. 
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Fig. 3-3. Z-component of solar force 
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Fig. 3-4. Magnitude of the solar radiation force 
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Fig. 3-5. Radial component of the solar radiation force 
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